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We discuss two topics in double parton scattering (DPS) theory that 
have been the subject of recent research interest. First, the role of 'double 
parton splitting' diagrams in DPS is discussed. We outline the 'dou- 
ble PDF' description of DPS, which was introduced a number of years 
ago. It is pointed out that under this framework, a certain structure is 
anticipated in the cross section expression for a 'double perturbative split- 
ting' diagram, which in the framework is regarded as DPS. We show that 
although this structure does indeed appear in the 'double perturbative 
splitting' diagrams, there is no natural reason to demarcate specifically 
this part of the graph as the DPS part, and indeed there is no natural part 
of these diagrams that can be regarded as DPS. There therefore appear to 
be some unsatisfactory features in the double PDF approach to describing 
DPS. The second issue we discuss is that of interference and correlated 
parton contributions to proton-proton DPS. We explain in simple terms 
why there can be such contributions to the DPS cross section. The poten- 
tial existence of such contributions was pointed out long ago by Mekhfi 
and more recently by Diehl and Schafer, but has been largely ignored in 
phenomenological investigations of DPS. 
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1 'Double Perturbative Splitting' Diagrams in 
Double Parton Scattering 

We define double parton scattering (DPS) as the process in which two pairs of par- 
tons participate in hard interactions in a single proton-proton (p-p) collision. DPS 
processes can constitute important backgrounds to Higgs and other interesting signals 
(see e.g. [1]), and can themselves be considered as interesting signal processes, since 
they reveal information about parton pair correlations in the proton. 

Making only the assumption that the hard processes A and B may be factorised, 
the cross section for proton-proton DPS may be written as follows: 

^iA,B) ^ X] y ]^'^a;a(i^6(Tifc^A(s = XiX3s)d-ji^Bis = X2XiS) (1) 
i,j,k,l a=l 

X Tij{xi, X2, b, Q% Q|)rH(a;3, x^, b; Q\, Q%) 

The cross section formula is somewhat similar to that used for single parton scat- 
tering (SPS), except that two parton-level cross sections a appear, and the PDF 
factors are two-parton generalised PDFs P (2pGPDs) rather than single PDFs. Note 
that in this formula the two 2pGPDs are integrated over a common parton pair trans- 
verse separation b - the transverse separation must clearly be identical in both protons 
in order that two pairs of partons meet in two separate hard interactions A and B. 
The DPS cross section cannot naturally be written in terms of PDFs individually 
integrated over their b arguments, as is the case for the SPS cross section. 

In many extant studies of DPS, it is assumed that the 2pGPD can be approxi- 
mately factorised into a product of a longitudinal piece and a (typically flavour and 
scale independent) transverse piece: 

Tij{x,,X2,b-Q\,Ql)c^D;^{xuX2;Q\,Ql)F{b) (2) 

Then, if one introduces the quantity (Je// via de// = l/[J F^b^d'^b], one finds 
that one may write crR^. entirely in terms of the longitudinal piece and (Jeff- 




In [2] a quantity D^J{xi,X2; Q^) is introduced, and an evolution equation for this 
quantity is given. We shall refer to the quantity and its evolution equation as the 
double PDF (dPDF) and the dDGLAP equation respectively. It is asserted in [3] 
that the dPDF is equal to the factorised longitudinal part of the 2pGPD in the case 
in which the two hard scales Q\ and Q\ are equal to a common value Q^. 

The dDGLAP equation contains two types of terms on the right hand side - 'in- 
dependent branching' terms corresponding to emission of partons from a pre-existing 
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Figure 1: (a) A diagram that apparently contributes to the leading order DPS cross 
section according to the framework of [3]. (b) The 'crossed box' graph. In this part of 
the figure, A and B are arbitrary single particle final states with Qi = Ql = > 0. 

pair, and 'single parton feed' terms corresponding to the perturbative generation of 
a pair from the splitting of a single parton. The single feed terms involve the leading 
twist single parton distributions as one might expect. Given this structure of the 
dDGLAP equation, with single feed terms included on the right hand side, a predic- 
tion of the framework suggested in [3] for calculating the proton-proton DPS cross 
section is that a part of the graph drawn in figure [U^a) should be included in the DPS 
cross section. The part that should be included is proportional to [log(Q^/A^)]"/(Te// 
at the cross section level, where A is some IR cutoff of order Aqcd, and n is equal 
to the total number of QCD branchings in figure [H^a) (including the two that only 
produce internal particles). This piece should be associated with the region of trans- 
verse momentum integration for the graph in which the transverse momenta of the 
branchings on either side of the 'hard processes' in the graph are strongly ordered. 

The question that then naturally arises is whether such a structure in fact exists in 
the cross section expression for the loop of figure [D^a). Starting from the conventional 
'Feynman rules' expression for the loop, it is not immediately obvious what the answer 
to this question is. Here we will focus on answering this question for the specific very 
simple 'crossed box' loop shown in figure [T]^b), which is predicted by the framework 
of [3J to contain a piece proportional to [log(Q^/A^)]^/(Te//. The issues raised in the 
treatment of this example carry over directly to the more general loop of figure [l](a). 

We expect the [log((5^/A^)]^/cre// piece in figure [Tt^b) to be predominantly con- 
tained in the portion of the cross section integration in which the external transverse 
momenta, as well as the transverse momenta and virtualities of the internal particles, 
are all small. This is actually the region around a certain Landau singularity in the 
loop integral known as the double parton scattering singularity In |5j, we obtained 
an analytic expression for the part of an arbitrary loop containing a DPS singular- 
ity associated with the loop particles emerging from the initial state particles being 
nearly on-shell and collinear, in the limit in which the external transverse momenta 
are small. Applied to the loop of figure Wljo) this reads (schematically, suppressing 
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helicity and colour labels): 



LDPSfigm^) j d^^^g^gqi^^ ^ - Q2)^g^qq{^ " X, -k) (4) 

X Mgq~,A{s = x(l - x)s)Mqg^B{s = x(l - x)s) + iq^q) 

In this formula, x = p2 ■ Qi/vi ■ P2, s = {pi + ^2)^, and k {Q2) is the component 
of k {Q2) transverse to the axis defined by the directions of the incoming particles. 
^g^qq{x, k) is the 0{as) light cone wavefunction (LCWF) to produce a qq pair from 
a [B] , with the quark having lightcone momentum fraction x and transverse momen- 
tum k with respect to the parent gluon. It can be factored into a A; and x dependent 
part, where the k dependent part is proportional to e • fc/fe^, e being the transverse 
part of the gluon polarisation vector. It is generally true that the k dependent part 
of the LCWF corresponding to any QCD splitting is proportional to 1/k. 

Inserting (jl]) into the standard 2 — )■ 2 cross section expression, and performing 
a number of changes of variable, we arrive at the following expression for the DPS 
singular part of the gg — ?■ AB cross section: 

o-DPS.figHIb) oc j Y\_dxidxiagg^Ais = XiXis)agg^Bis = X2X2S) 



(5) 



1=1 



X / 7:r^'^g^qq{xi^X2,r)Tg^gg{xi,X2,-r) 



(2vr)^ 

as ■■ rk'<oiQ'') [^^i^],[^_i^p- 

Tg^gg{xi,X2,r) oc—S{l-Xi-X2)T^{xi,X2) d k | | (6) 

ATT J [k+ [k 

T''^{xi,X2) contains a function of Xi and X2 that may be regarded as a '1 — )■ 2' 
splitting function, multiplied by a constant matrix in transverse spac^. r is equal 
to the transverse momentum imbalance of one of the quarks/antiquarks in the loop 
between amplitude and conjugate, and is the Fourier transform of the parton pair 
separation b in the qq pair emerging from either gluon. Tg^gglxi, X2, r) can therefore 
be thought of as the 0{as) transverse momentum-space 2pGPD to find a qq pair 
inside a gluon. Note that the expression here effectively coincides with that of in 
which a cross section expression for the box of [T](b) is obtained starting from a pure 
DPS view of the box. 

Let us consider the part of the integral ([5]) that is associated with the magnitude 
of the imbalance r being smaller than some small cut-off A that is of the order of 



*Note that the cross section is reahy a sum of terms with different T'^^ {xi,X2) factors in the 
g qq 2pGPDs. This is associated with the fact that, from the point of view of the quarks, there is 
an unpolarised diagonal contribution to the process plus polarised and interference contributions in 
colour, spin, and flavour space. See section [2] for a discussion of correlation and interference effects 
in DPS processes. 
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^QCD- The contribution to the cross section from this portion contains a log^(Q^/A^) 
factor multiphed by (which can be thought of as an effective 'l/cTg//' factor for 
this contribution). The majority of this contribution comes from the region in which 
the transverse momenta and virtuahties of the quarks and antiquarks in the gg — )■ 
AB loop are much smaller in magnitude than \/Q^ (i.e. the region in which the 
assumptions used to derive (jl]) apply), which is a necessary feature of a contribution 
to be able to regard it as a DPS-type contribution. By making a specific choice of A 
(let us call this A^), one could obtain an expression which is exactly in accord with 
the expectations of [3] - that is, a product of two large DGLAP logarithms multiplied 
by the same l/o"e// factor that appears in diagrams in which the parton pair from 
neither proton has arisen as a result of one parton perturbatively splitting into two. 
The l/o"e// factor for these diagrams presumably has a natural value of the order of 
1/-Rp that is set by the nonperturbative dynamics {Rp = proton radius). 

The fact that we have to make a somewhat arbitrary choice for A in order to arrive 
at the result anticipated by the framework of [3j is concerning. There is nothing in 
the calculation of figure lH^b) to indicate that we should take the region of it with 
|r| < As as the 'DPS part' - the scale A^ does not naturally appear at any stage 
of the calculation. There is no more justification for taking the part of the box with 
\r\ < As to be the DPS part than there is for, say, taking the piece with \r\ < 2As, 
or that with |r| < As/2, to be the DPS part. 

There therefore appear to be some unsatisfactory features of the framework of [3] 
with regards to its treatment of the box in figure [U^b). It is clear that these issues 
will also be encountered for the case of the arbitrary 'double perturbative splitting' 
graph in figure [U^a). One obtains a result that is consistent with the framework of [3] 
if one demarcates the portion of the cross section integral in which the transverse 
loop momentum imbalance between amplitude and conjugate is less than As as DPS, 
but there is no natural reason to do this. There is no distinct piece of figure [U^a) 
that contains a natural scale of order Aqqd and is associated with the transverse 
momenta inside the loop being strongly ordered on either side of the diagram - most 
of the contribution to the cross section expression for this graph comes from the 
region of integration in which the transverse momenta of particles inside the loop 
are of 0(\/Q^) (except perhaps at low x values for the 'hard subprocesses' in the 
graph [8j). It is therefore perhaps the case that we should not regard any of this 
graph as DPS. Treating the graph in this way has the advantage that we do not 
perform any double counting between DPS and SPS - the graph of figure [It^a) is in 
principle also included in the SPS pp — )■ AB cross section (albeit as a very high order 
correction that will not be included in practical low order calculations, if the number 
of QCD emissions from inside the loop of the graph is large). 

One can gain some insight into the source of the problems in the framework of [3] 
by looking at the 6-space 2pGPD corresponding to ([6]). This comes out as being 
proportional to 1/6^ - this behaviour (which was first spotted in |7j) can be traced 
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to the fact that the g ^ qq LCWF in b space (hke any LCWF corresponding to a 
QCD perturbative sphtting) is proportional to 1/6, and r(6) ~ Note that this 

behaviour is very different from the behaviour of all 2pGPDs that is anticipated by 
the dPDF framework (i.e. smooth function of size Rp). There is no natural feature in 
the product of two 'perturbative splitting' 2pGPDs that is of size Rp and can be nat- 
urally identified as DPS. A key error then in the formulation of the dPDF framework 
seems to be the assumption that all 2pGPDs can be approximately factorised into 
dPDFs and smooth transverse functions of size Rp. A sound theoretical framework 
for describing proton-proton DPS needs to carefully take account of the different 
b dependence of pairs of partons emerging from perturbative splittings, whilst si- 
multaneously avoiding double counting between SPS and DPS. For recent proposals 
from other authors as to how proton-proton DPS should be described theoretically, 
including further discussion of the 'double parton splitting graphs', see [8l|9]|ll 

2 Interference and Correlation Effects in DPS 

It was pointed out long ago by Mekhfi [TT] that there can be contributions to the 
p-p DPS cross section associated with spin and colour correlations and interference 
effects in spin and colour space, even when the colliding protons are unpolarised. The 
issue was taken up again recently by Diehl and Schafer [7] who demonstrated that the 
correlation and interference contributions may be sizeable, and made the observation 
that interference effects in flavour space can also contribute to unpolarised p-p DPS. 
In this section we explain in simple terms why there can be interference and correlated 
parton contributions to the unpolarised p-p DPS cross section, where there are no such 
contributions to the corresponding SPS cross section. We hope that this explanation 
may be of aid to those less familiar with the subject, and refer the reader to jTlflT] 
for more details. 

We recall that the cross section for leading twist single parton scattering processes 
is calculated from 'cut diagrams' with the structure of figure |2I^a). For definiteness we 
have taken the SPS process to be Drell-Yan in the figure, but the details of the final 
state are not important for our discussion. Now, if we consider the parton 'returning' 
to (say) the bottom proton on the right hand side of the diagram, then we see that 
it must have exactly the same flavour and colour as it 'left' with on the left hand 
side. This must be the case otherwise it cannot 'reform' the original proton when 
it combines with the spectators on the right hand side. So there can be no flavour 

^It is worth pointing out in passing that there is perhaps a double scattering process which does 
directly involve the dPDF of the proton. This is the contribution to proton-heavy nucleus DPS 
associated with partons from two separate nucleons interacting with two partons from the proton. 
The reason why this probes the dPDF, and p-p DPS does not, is that in this case the 'probe' parton 
pair coming from the nucleus has a (roughly) flat distribution in b, whereas in p-p DPS some parton 
pairs (i.e. perturbatively generated ones) do not. For more details and a discussion of how the 
two-nucleon contribution to proton- heavy nucleus DPS might be extracted experimentally, see [lOj . 
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Figure 2: (a) Leading twist diagram for single Drell-Yan (as an example process) in 
proton-proton collisions, (b) A diagonal in colour (left) and colour interference (right) 
diagram contributing to the double Drell-Yan DPS process. 

and colour interference contributions to p-p SPS. When the colliding protons are 
unpolarised, symmetry forbids any contribution to the SPS cross section associated 
with helicity or transversity polarisation effects. For similar reasons, there cannot 
be any contribution to the SPS cross section associated with the analogous effects in 
colour space. The only PDFs that contribute to the unpolarised SPS cross section 
are therefore the unpolarised diagonal colour-summed PDFs. 

The cross section for DPS processes is calculated from cut diagrams with the 
structure of figure E^b) in which two partons 'leave' each proton on the left, inter- 
act, and then 'return' on the right. In this case, the fact that the proton must be 
reformed at the end only imposes constraints on the 'sum' of the two partons' quan- 
tum numbers. For any DPS process, the possibility arises for interference diagrams 
to contribute in which one or more of the discrete quantum numbers get swapped 
between the partons before they return to the proton - provided that a swap in the 
opposite direction happens in the other proton (an example of a colour interference 
diagram that contributes to double Drell-Yan is the right hand graph of figure |2]^b)). 

There are also contributions to the DPS cross section associated with polarised 
2pGPDs that can be nonzero even when the colliding protons are unpolarised. The 
reason for this is that, with two partons participating from each proton in DPS, 
there can be effects relating to the correlations in spin between the active partons 
themselves. So, if for example the chance to find a pair of quarks in the proton 
with their helicities aligned differs from that to find a pair of quarks with opposing 
hehcities, then AgiAg2 = Q'lt Q'2t Q2i —Qii Q'2t ~Q'it I2i7^ 0. In a similar way 
there will be contributions to p-p DPS associated with colour correlations. 

Despite it being pointed out long ago that p-p DPS may be affected by interference 
and correlated parton effects, such effects are generally not considered in phenomeno- 
logical analyses of this process. More detailed studies need to be performed on these 
effects, including an examination of the effect of evolution on them (to what extent 
does evolution 'wash out' correlations?) and it would be desirable to have more re- 
fined estimates of the size of interference and correlated parton contributions than 
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were performed in [7]. One would require low-scale inputs for the interference and 
correlated parton two-parton distributions in order to make such estimates (along 
with the requisite evolution framework), very approximate forms for which could per- 
haps be extracted from proton models. Such forms would of course not be reliable 
at low X owing to the fact that one cannot fit parton densities at low x even at low 
without including a number of 'nonperturbative' gluons and sea quarks [12], and 
proton models typically only include the lowest few Fock states. An alternative ap- 
proach for obtaining 'first guess' inputs for some of the distributions via single-parton 
GPDs is given in [7]. 



References 



[1] 

[2] 

[3] 
[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 

[11] 
[121 



A. Del Fabbro and D. Treleani, \Phys. Rev. D61 (2000) 077502, 
[arXiv : hep-ph/99 1 1358} 

G. M. Zinovev, A. M. Snigirev, and V. P. Shelest, 



\Theor. Math. Phys. 51 (1982) 523-528 



A. M. Snigirev, Phys. Rev. D68 (2003) 114012, [arXiv:hep-ph/0304172 



Z. Nagy and D. E. Soper, 'Phys. Rev. D74 (2006) 093006 
,arXiv:hep-ph/0610028^ [hep-ph] . 



J. R. Gaunt and W. J. Stirling, \JHEP 1106 (2011) 048 
.arXiv: 1103 .1888 fhep-phl , 



A. Harindranath, R. Kundu, and W.-M. Zhang, \Phys.Rev. D59 (1999) 094013 
|arXiv:hep-ph/9806221, [hep-ph], 



M. Diehl and A. Schafer, Phys . Lett. B 698 (2 011) 389-402 
.arXiv: 1102.3081. fhep-phl. 

M. Ryskin and A. Snigirev, \Phys.R ev. D83 (2011) 114047| 
larXiv: 1103.3495. rhep-pFTl 

B. Blok, Y. Dokshitzer, L. Frankfurt, and M. Strikman, 
larXiv: 1106.5533. fhep-phl [ 

M. Strikman, D. Treleani, Phys. Rev. Lett. 88 (2002) 031801. 
|arXiv:hep-ph/0111468| [hep-ph]] . 

M. Mekhfi, \ Phys.Rev. D32 (1985) 2380 



M. Gluck, P. Jimenez-Delgado, E. Reya, Eur. Phys. J. C53 (2008) 355-366. 
[arXiv:0709.0614l [hep-ph] ] . 



7 



